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1 Introduction 

The name spiked oscillator generally refers to Hamiltonians with potentials 
containing a parabolic term x 2 describing ordinary quantum oscillator, and 
a singular term l/x a forming a sharp spike at x = 0. A typical Hamiltonian 
of a spiked oscillator is of the form given by Harrell||: 

d o A 
H = — — + x H < x < oo. 

dx 1 x a 

It has found variety of uses in atomic, molecular, nuclear and particle physics 
since it provides the simplest realistic model of interaction potentials due to 
its repulsive core x~ a . 

Some forms of spiked oscillators have exact solutions, such as the one 
described by one dimensional Hamiltonian, examined four decades ago by 
Goldman and Krivchenkov[|J 

d z x a 

whose exact solutions 4> n {x) £ [0, oo) satisfy Dirichlet's boundary condition 

0n(O)=O. 

Models with exact solutions can serve as starting points for analyses of 
related cases where exact solutions do not exist and where a perturbative 
analysis is required A good example can be found in the papers of Hall 
et al[[ij]p|], who slightly modified the above Hamiltonian, computed a basis 
given by a set of exact orthonormal eigenfunctions (f) n (x) for L 2 (0, oo), and 
then used this basis in variational analysis of the Hamiltonian 

d n 2 A A 

They argued that their singular basis forms much better starting point for 
perturbation techniques than the basis of the ordinary harmonic oscillator, as 
used in earlier work of Aguilera-Navarro et al[||, who had applied a perturba- 
tive variational analysis to the lowest eigenvalue of the Harrel's Hamiltonian. 

In this context, we present an analysis of family of following spiked Hamil- 
tonians 
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parametrized by discrete values of I = 0,1,2... Depending on interpreta- 
tion of variable x, this form can either describe a family of one dimensional 
linear oscillators, or a family of N dimensional radial oscillators for odd 
N = 1,3,5.. What is traditionally known as the three dimensional isotropic 
radial oscillator, with known exact solutions, can be reduced to a special 
case of N — 1. We shall show that solutions for radial oscillators in even 
dimensions N = 2,4,6.. do not form well defined orthogonal bases and, as 
such, they are of no practical interest. 

The solutions we provide are exact and much simpler than any of those 
described in this introduction, and therefore they are good candidates for 
orthonormal bases that could be used in perturbative analyses of oscillators 
with spiked potentials other than x~ 2 . 

The factorization formalism developed here is completely algebraic and 
somehow related to methods of Superstring Quantum Mechanics (SUSY 
QM). Modern algebraic factorization methods are overviewed in the paper of 
Rosas-Ortiz and closely related SUSY methods are described for example 
in introductory parts of the papers by Junker and Roy fl I ] [12], and Levai 
and Royp~3"l . 



The common premise of all those methods is the use of algebraic recursive 
manipulations on two Hamiltonians: one with known spectrum and known 
eigenvectors and another which is isospectral to the former but whose eigen- 
vectors are to be found. The methodology developed here uses the same 
premise, but goes a bit further by manipulating not just two Hamiltonians 
but the uncountable set of related Hamiltonians (fl|). On the other hand, 
the SUSY and similar methods rely on general solution of underlying Ri- 
catti equation, which is completely ignored in the approach presented in this 
paper. 

Quantum models with exact solutions play an important role and there 
is a renewed interest in enrichment of the traditional narrow set of exactly 
solvable models, such as harmonic oscillator, hydrogen-like potential, Morse 
potential and square well potential. It is only fair to mention other re- 
lated methods, which admit certain type of non-hermitian Hamiltonians with 
complex-valued potentials giving rise to real energy spectrum. For example, 
Cannata et al[|10|| use Darboux method to derive a general class of complex 
potentials with energy spectra identical to that of regular harmonic oscilla- 
tor. Bender et al|||| replaces the requirement of Hamiltonian hermiticity 
by weaker requirement of PT- symmetry , which is satisfactory condition for 
real energy spectra. They examine complex harmonic and anharmonic os- 
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cillators, complex square wells, etc. This avenue is being explored by many 
other researchers. 

To complete this introduction, we just mention in passing that exactly 
solvable problems can be categorized in one of three classes: Exactly solvable 
(ES), Quasi exactly solvable (QES) and Conditionally exactly solvable (CES). 

In the next section we present an overview of our findings, while the 
remaining sections provide specific details. 

2 Outline 

The first Hamiltonian, Hq, from the family of discretely spiked oscillators ([!]) 
describes the well known harmonic oscillator, without the spike, of known 
discrete energy spectrum and known ladder operators a and aj. Aside from 
this simple case with parabolic potential, the potentials of remaining mem- 
bers of this family have strong discontinuities at x = and do not resemble 
the first case at all. Yet they all are isospectral; that is, they all have the dis- 
crete energy spectra identical to that of the simple harmonic oscillator. The 
only difference is in the definition of their ground state energies, Eiq = 21 + 1, 
in units of \ujTi. 

We have chosen to denotate the coresponding ground states as |0,0), 
|1, 0), |2, 0), etc. for I = 0, 1, 2.., respectively. Table |] visualizes the relation- 
ship between energy states and energies for the entire family of discretely 
spiked oscillators. 

For each Hamiltonian Hi a pair of intertwining operators can be defined 

bj +1 \l, k) oc |Z + l,fc) 
& {+ i|Z + l,fc) cx \l,k) 

which connect eigenstates of two 'neighbours', Hi and Hi + \. Both represent 
actions coresponding to diagonal traversal of the Table ^ in NE and SW 
directions, respectively. 

In addition to intertwining operators, we can also define the ladder oper- 
ators ai and a\ that allow for traversing of the Table [I] in vertical direction. 
The annihilation operator changes the state \l,k) with energy Ei^ to the 
state \l, k — 1) of lower energy Ei >k _i 

ai\l,0) = 

a t \l,k) = ai\l,k — l). 

Similarly, the creation operator changes the state \l,k) with energy Ei t k to 
the state \l, k + 1) of higher energy E^ k+1 

a\\l,k) = a\\l,k + 1) 
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Table 1: States and energies of discretely spiked harmonic oscillators. Some 
states are unphysical and will be later rejected. 



1 = 


I = 1 


I = 2 


I =... 


Energy 










\2,k) 






2/b + 5 






\hk) 




\2,k- 


1) 




2fc + 3 


\0,k) 




\l,k- 


1) 


2, k - 


2) 




2fc + 1 


jo, A: — 


1) 


\l,k- 


2) 








2A; - 1 


jo, A; — 


2) 












2k -3 










|2,2) 






9 






|1,2) 




|2,1) 






7 


|0,2) 




IM> 




|2,0) 






5 


|o,i) 




|i, o> 










3 


|0,0) 














1 



where a\ and a[ are the normalization constants. 

Known ladder operators do and a\ for the classical harmonic oscillator are 
used to recursively prove the existence of similar operators for the remaining 
Hamiltonians Hi. The are given by formulas 

a i+i = bj +1 aibi + i 

The condition ai\l, 0) = allows to generate explicit form of the state 
|/,0) in position representation; that is the vaweform function 4>i t o(x) of the 
ground state. For I = 0, this is the bell-shaped function, o ,o oc exp (^f-), 
subject to normalization. Having computed the ground waveform 0/ jO) one 
can use the creation operator a\ to recusively generate all other eigenstates 
\l,k), or coresponding waveforms (j>i^{x). 

This seemingly completes the scheme since, thanks to the pairs of ladder 
and intertwining operators, the entire Table [I] can be created - starting from 
the left lower corner and zigzagging one step a time, with or without the help 
of ladder operators a\, a' 2 , etc. 

However, we must also require that each set of state vectors is a basis in 
its own space; that is, they all are mutually orthogonal and normalized 

(l,j\l,k) = 5 jk . 
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But some of the wavefunctions are singular at x — 0, since they contain 
factors -, -2, etc. Depending on whether the Hamiltonians ([I]) describe 
linear or radial oscillators and, in the latter case, depending on the space 
dimension, some such functions are not square-integrable and therefore they 
must be rejected as unphysical. Those wave functions that are physical must 
additionally pass the test of orthogonality in order to qualify as members of 
useful bases. To take it all into account we therefore specialize the generic 
scheme and draw additional conclusions specific to dimensionality of the 
problem. All of this is supported by a Haskell program, which tests the 
theory and provides useful and flexible computational tool. 

In the following sections we shall present a detailed account of this outline. 

3 Operators 

We shall define two pairs of state generation operators: intertwining oper- 
ators, affecting the quantum number I, and ladder operators affecting the 
qunatum number k. The former correspond to diagonal traversal of Table p] 
and the latter - for the vertical traversal. 

3.1 Intertwining operators 

Let's define two operators 6/ and its adjoint b\ 



d 

dx 



+ AO), 



(2) 




where 



(3) 




(4) 




and, as a consequence, the commutator [bi,b\] is not zero either, since 



dx 2 
d 2 



+ /3f(x)-i3' l (x) = H l + (2l-l), 



(5) 




dx 2 



+ A 2 (x)+A(x)=^_ 1 + (2/+l), 



(6) 
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where 

d 2 

H_ x = # = -— + x\ 

dx 1 

and where the remaining Hamiltonians are defined by (|IJ). 

With the help of the above equations we can establish two very useful 
formulas, which we will use later: 

Hfij-blH^ = 2b\, (7) 
Hi-ibi-bM = -2b h (8) 

These are quasi commutative relations, connecting two "neighbouring" 
Hamiltonians. Let us now consider specific cases of /3i(x). 

3.1.1 Ordinary oscillator 

The case I = is very special since 

b% = # -l, 
b bl = #_x + l 
= #0 + 1- 

Both equations refer to the same Hamiltonian of ordinary harmonic oscil- 
lator, and the operators b and 6 degenerate to the familiar ladder operators 



a 



d 



bl = —r + x, (9) 



dx 



a = b = — + x, (10) 
dx 



which have these important properties: 



[oqA] = 2, (11) 
[# ,a ] = -2a , (12) 
[# ,a ] = 2al (13) 

Consequently, the energies are given by Eq^ = 2k + 1 for k = 0, 1, 2..., in 
units of kuih. We will skip other details, since this is the very known case. 
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3.1.2 Discretely spiked oscillators 

The interesting part begins here. First, let us prove that b] +1 \l,k) is an 
eigenvector of Hamiltonian Hi + \ by assuming that this eigen equation holds 

H^l.k) = E Uk \l,k) 

and making use of formula (|7|): 

H l+1 \bj +1 (l,k)) = H l+1 b} +1 \l,k) 

= b\ +1 H l \l 1 k) + 2b\ +1 \l,k) 
= (E l>k + 2)\b} +1 (l,k)). 

The operator b\ +1 transforms the state \l,k) with energy E^ k into the state 
\l + l,k) with energy E l+ljk = E l)k + 2 

\l + l ,k) = — r 1 bj +1 \l,k). (14) 
1 1 V2k + 4/ + 4 1+1 1 ' V ' 

This coresponds to the diagonal move in NE direction in Table [I]. One can 
verify the normalization factor in equation (]14|) by comparing the norm of its 
both sides, making use of (^|) and noticing that the energy of 1-th oscillator 
in the state \l,k) is given by 

E ltk = 2(1 + k) + l. (15) 

The latter can be easily established with the help of Table [I]. 

Similarly, we can prove that state bi + i\l + 1, k) is the eigenvector of Hamil- 
tonian Hf. 

Hi{b l+1 \l + l,k)) = {b l+1 H l+1 -2b l+1 )\l + l,k}) 
= (E l+1>k -2)(b l+1 \l + l,k)). 

The operator b i+ i transforms the state \l + l,k) with energy Ei + i )k into 
the state \l, k) with energy Ei k = E [+ i )k — 2 

\l,k) = -;= bi +1 \l + l,k). (16) 

1 ' V2k + 4Z + 4 + 1 1 y 1 

This coresponds to the diagonal move in SW direction in Table 0. 
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3.2 Ladder operators 

We shall now define two operators 

4+i = h \+A h i+i ( 17 ) 

and 

a i+ i = b] +1 aik + i (18) 

and prove that they indeed have the expected properties of ladder operators. 
We shall first prove by induction, that 

[H h a\]=2a\ (19) 

This holds true for the case I = 0, as seen in (13). Assuming now, that the 
above is true for I — 1, the proof is as follows: 

E x a\ = HfijaUk 

ij from 
= (bjH l _ 1 + 2bl)al 1 b l 

f from Q 
= b}(a}_ 1 H l -. 1 + 2a}_ 1 )b l + 2a} 

•ft" from by induction 
= bjaj^H^bi + 4aJ 
= b\a\_ l {b l H l -2b l )+Aa\ 

f from |] 
= bja\_ 1 b l H l + 2aj 

= 4 H i + 2 4 
ij from 

Similarly, we can prove by induction that the following holds 

[H h ai] = -2ai. (20) 

It immediately follows from (|19"D and ( p0[) that operators 4 an d a i create 
and annihilate eigenstates of Hamiltonian Hi, respectively 

Hi4\l,k) = a\Hi\l,k)+2al\l,k) 
= (E ljk + 2)4\l,k), 

H t ai\l,k) = a t Hi\l,k) -2a t \l,k) 
= (E ltk -2) ai \l,k). 
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Finally, we can prove by induction that if do 1 0,0) = then similar relation 
holds for any other oscillator (/ > 0): 

a,|Z,0> = 0. 



4 Wave functions 

4.1 Haskell module Spike 

The remaining part of this paper refers to results obtained with the help of 
a little program, Spike.hsj^ written in functional lazy language HaskellQ. The 
program is used for generation of wavef unctions, for normalization of solu- 
tions in n-dimensional spaces and, generally, for testing the theory presented 
in the previous sections. The program performs some algebraic manipula- 
tion of Laurent series, such as computation of derivatives, integrals, sums 
and products, and directly implements raising and intertwining operators 
previously discussed. It is best used in interpreting environment^, since one 
can evaluate any formula from module Spike in any order of one's choice. 
Module Spike relies on other standard Haskell modules and on our module 
Fraction. hs PI. For interesting introduction to Haskell as a coding tool for 



scientific applications see|L4 . 



4.2 Laurent series 

Any analytic function f(z) in an open annulus r < \z — Zq\ < R can be 
expressed as the sum of two series 

oo oo 
3=0 j=l 

Such an expansion, containing negative as well as positive powers is called 
Laurent series for f(z) in this annulus. In Haskell we will represent it as a 
pair of two lists 

f = L [a , a 2 , •••] [a_i, a_ 2 , ...], 

where coefficients of expansions can be any numbers belonging to class Num, 
such as Integer, Double, Complex Double, etc. 



lr The program is available at [http : / / www. numeric-quest .com /haskcll / Spike . hs 
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See home page of Haskell at http://www.haskell.org 



3 Haskell interpreter Hugs is available for free at http:/ /www. haskell.org/hugs/. It runs 
on any major platform. 



Module Fraction is avauilable at http: / /www. numcric-qucst.com/haskcll/Fraction.h4 
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The following defines the type of data structure representing Laurent 
series 

data Laurent a = L [a] [a] , 

where a is a type variable of some numeric kind. It is convenient to choose it 
as Integer (of unlimited size) since this leads to highly accurate results. For 
computations requiring normalization the data Fraction will be used instead 
of traditional Double since this guarantees very good accuracy, with relatively 
small performance loss. 

As it will soon become clear, wavefunctions of spiked oscillators can be 
represented as products of two functions 



where fi,k{x) is a real function of x and can be expanded as Laurent series 
of type Laurent Integer . 

4.3 One step operators 



results in manipulation of two Laurent series pi and fk,i, where pi(x) is a 
short Laurent series, represented in Haskell as L[0, 2] [I]. The expression 
{—fki{x) + Pi(x)fk,i{ x ) ^ s implemented in Haskell as function b'. It takes 
as an input an integer I and a Laurent series fi^ of type Laurent Integer , 
and then produces another Laurent series by invoking few primitives, such 
as multiplication, addition and differentiation defined for such series. Notice 
the gcd normalization of the results, which keeps the size of integers in check 
and also provides a canonical standard for comparison of different methods 
of computations. 

b ; : : Integer -> Laurent Integer -> Laurent Integer 

b' 1 f = gcdNormalize (p * f - diff f) 




(21) 




where 



p = L [0,2] [1] 
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Implementation V I of operator bj is the only Haskell function needed to gen- 
erate all eigenfunctions for all spiked harmonic oscillators. But to comple- 
ment the theory developed in the previous sections, below are the remaining 
operators. Firstly here is function bl coresponding to intertwining operator 
k 

b : : Integer -> Laurent Integer -> Laurent Integer 

b 1 f = gcdNormalize (p * f + diff f) 
where 

p = L [] [1] . 

Next is the implementation of the rasing operator a): 

a' : : Integer -> Laurent Integer -> Laurent Integer 
a' = b' 

a' 1 = b' 1 . a' (1-1) . b 1 

Notice that function a' = b' is a special case. It coresponds to the raising 
operator a = b\ - exactly as specified in the theoretical sections. Finally 
here is the implementation of the lowering ladder operator af. 

a : : Integer -> Laurent Integer -> Laurent Integer 
a = b 

a 1 = b' 1 . a (1 - 1) . b 1. 
4.4 Cumulative operators 

Since the one-step operators from the previous section are to be applied re- 
cursively to some base functions it is convenient to define cumulative versions 
of these operators. 

The cumulative intertwining operator, twine' I, is a composition of (NE) 
operators bi] &;_if ... &if 

twine' : : Integer -> Laurent Integer -> Laurent Integer 
twine' = id 
twine' 1 = b' 1 

twine' 1 = b' 1 . twine' (1 - 1) 

The cumulative raising ladder operator, ladder' I k, is a composition of k 
raising ladder operators a\ a\ a\... 

ladder' : : Integer -> Integer -> Laurent Integer 

-> Laurent Integer 
ladder' 1 = id 
ladder' 1 1 = a' 1 

ladder' 1 k = a' 1 . ladder' 1 (k - 1) . 
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4.5 Generating wavefunctions 

Using these tools we can easily generate sets of eigenfunctions for spiked 
oscillators. When / = the first wavefunction is given by e~ x ' 2 and that 
implies that / ,o is simply 1, and its corresponding Laurent series is L [1] []. 
Here [] represents an empty list - meaning no negative powers at all. By 
executing 

iterate (a' 0) (L [1] [] ) 

we can generate infinite list of gcd-normalized eigenfunctions for ordinary 
harmonic oscillator. If we want to retain six, say, wave functions we need to 
cut this infinite list down to six elements 



take 6 (iterate W 0) (L [1] [])). 



A slightly prettiefied output of the above looks like this: 

— Meaning: 

L [1] [], — 1 

L [0,1] □ , — x 

L [-1,0,2] [] , 1 + 2x~2 

L [0,-3,0,2] [] , 3x + 2x~3 

L [3,0,-12,0,4] [], — 3 - 12x~2 + 4x~4 

L [0,15,0,-20,0,4] [] — 15x - 20x~3 + 4x~5. 



] 



It is evident that these solutions do not contain negative powers of x. In fact 
they represent (gcd normalized) Hermite polynomials in disguise, and this is 
exactly what one would expect from ordinary harmonic oscillator. To test 
implementation of the lowering ladder operator a , execute 

take 6 (iterate (a 0) (L [0,15,0,-20,0,4])). 

This should recreate the above list in reverse. In addition, the expression 

a (L [1] []) 

generates a special kind of Laurent series 
L [] [], 

which is the Laurent series for function zero. 

To find eigenfunctions of the first spiked oscillator (/ = 1) we have two 
choices. Firstly, we can apply operator b\, or function b' 1 to every eigenfunc- 
tion of ordinary oscillator, by executing 
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map (b> 1) (take 6 $ iterate (a' 0) (L [1] [])). 

Alternatively, we can iteratively apply function a' 1 to the first eigenfunction 
/i,o — bl(L [1] []) of the oscillator I = 1 

take 6 $ iterate (a' 1) f where f = b' 1 (L [1] [] ) . 



Either way the output is as follows 









— Meaning: 






L 


[0,2] 


[1], 


— 2x + 1/x 




L 


[0,0,1] 


□ , 


— x~2 






L 


[0,-4,0,4] 


[-1], 


4x + 


4x~3 


- 1/x 


L 


[0,0,-5,0,2] 




5x~2 + 


2x~4 




L 


[0,18,0,-36,0,8] 


[3], 


— 18x 


36x"3 


+ 8x"5 


L 


[0,0,35,0,-28,0,4] 


[] 


— 35x~2 - 


28x~4 


+ 4x~6 



] 



This time, every other solution contains term with 1/x, and as such it might 
or might not be acceptable as physical solution, depending on the space 
dimension under consideration. More about it in later sections. 

To generate similar list for the second spiked oscillator (I = 2) we need 
to apply two functions b' 1 and b' 2 in succession to a list of eigenf unctions of 
ordinary harmonic oscillator 

map (b> 2 . b' 1) (take 6 $ iterate (a' 0) (L [1] [])), 

or, alternatively, apply the raising ladder function a' 2 to the lowest eigen- 
function of the oscillator / = 2 

take 6 $ iterate (a' 2) f where f = (b> 2 . b' 1) (L [1] [] ) . 

Both methods lead to the same output: 
[ 



L 


[4,0,4] 


[0,3] , 


L 


[0,0,0,1] 


[], 


L 


[-6,0,-12,0,8] 


[0,-3] , 


L 


[0,0,0,-7,0,2] 




L 


[24,0,72,0,-96,0,16] 


[0,9] , 


L 


[0,0,0,63,0,-36,0,4] 


[] 



] . 



Here is a general formula for generation of infinite lists of wavefunctions for 
spiked harmonic oscillators I = 0, 1, 2.. 
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wavef unctions : : Integer -> Laurent Integer 

wavef unctions 1 = map (twine' 1) $ iterate (a' 0) (L [1] [] ) 

To generate just one eigenfunction fi :k the following Haskell function can be 
used 

wavefunction : : Integer -> Integer -> Laurent Integer 

wavef unction 1 k = twine' 1 (ladder' k (L [1] [])) 

5 Normalization and orthogonalization 

We have shown that the solutions for spiked oscillators are of the form 0^ = 
fi,k( x ) e ~ x I 2 1 where fi,k( x ) are rea l functions of x and, specificly, they are 
some Laurent series - admitting both positive and negative powers of x. For 
a specific case of ordinary oscillator, these are the Hermite polynomials, with 
positive powers only. 

The scheme presented so far accepts all solutions, but does not check 
whether they are square-integrable and, in addition, whether they obey the 
orthogonality relations 

= Sij. 

Due to presence of the term - in the definition of b\ some of the wavefunctions 
(fii^k are tainted by powers of - and as such might not be admissible as physical 
solutions for a dimension under consideration. Terms such as 

x 

when considered in one dimension, imply normalization integrals such as 

f + °° 1 _ 2 

/ ~2 e ~ X dX 
J — oo X 

which are not integrable. On the other hand, if variable x is treated as three 
dimensional radius then the normalization integral like this 

r+oo 1 2 9 

/ ~^e~ x 4irx dx, 
Jo x 2 

is integrable. 
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5.1 Admissible physical solutions 

A general pattern of physically acceptable solutions is as follows: 



• All odd eigenfunctions 4>i,k for A; = 1,3, 5... of any oscillator / = 0, 1, 2... 
are physically acceptable, regardless the space dimension. 

• All even eigenfunctions <j)^k for k — 0, 2, 4... of any oscillator I — 0, 1, 2... 
are acceptable as physical solutions only when N > 21 + 1, where N is 
the space dimension. 

The above statements can be verified by Haskell function physicalPattern, 
which produces infinite list of booleans - each stating whether a solution k 
for oscillator I and space dimension N is admissible as physical solution 

physicalPattern : : Integer -> Integer -> [Bool] 
physicalPattern = map (isPhysical n) (wavef unctions 1) 

isPhysical : : Integer -> Laurent Integer -> Bool 
isPhysical n f 

I length us' == = True 

I otherwise = False 

where 

L us us ' = f * f * dV 

dV = diff $ L ((take (f romlntegral n) 

$ repeat 0)++[l]) [] 

For example, even solutions for oscillator / = 2 are unphysical in three 
dimensional space, N = 3, as seen below (Haskell lists are traditionally 
indexed from zero) 

take 6 $ physicalPattern 3 2 

==> [False, True, False, True, False, True] . 

5.2 Linear normalization 

In order to normalize the eigenfunctions of one dimensional linear spiked os- 
cillators and to test their orthogonality conditions we need algebraic formulae 
for the values of improper integrals 



Pr 



n 




dx. 



They are as follows 
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Pi 







n 




n-2- 



Notice that integrals with odd powers of x n are all zero. 

All solutions for the ordinary harmonic oscillator (I = 0) are physically 
acceptable since their Laurent series do not contain coefficients of negative 
powers of x. What's more, they form the orthogonal system of eigenf unctions 



Physically acceptable solutions for the remaining, spiked, oscillators (I > 
0) are those with odd values oik, k — 1, 3, 5... Each oscillator I has its own 
orthogonal basis. This can be easily checked with the help of the Haskell 
function bracket, the scalar product of two (not normalized) wavefunctions - 
each annotated by a pair of integers (I, k). 

bracket :: Integer -> (Integer , Integer) 

-> (Integer, Integer) -> Fraction 
bracket n (1, k) (1' , k') 

= scalarProduct n f g 

where 

f = wavef unction 1 k 
g = wavef unction 1' k' 

The integer n specifies the space dimension n = 1,2,3.. when computing 
volume integrals for radial cases. But we can use the same function for 
computing linear integrals as well - by setting, somewhat artificially, this 
argument to zero, n — 0. For example 

decimal 8 $ bracket (7,1) (7,1) ==> 14034.40729347 
decimal 8 $ bracket (7,1) (7,5) ==> 

Notice the convertion from fractional to decimal representation, specified 
here with accuracy to eight decimal places. All internal computations are 
performed on integers of unlimited size and on fractions made of such inte- 
gers. The accuracy can be extremely good, depending on the value of eps, 
set at the top of module Spike. It is only for presentation purposes that we 
convert fractions to decimal representation. 

The following table summarizes the case of linear spiked oscillators 



(0,k\0,k')=8 w 
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/ = 





k = 


0,1,2... 




= 1,3,5... 


I = 


1 


k = 


1,3,5... 




= 5,9,13... 


I = 


2 


k = 


1,3,5... 


E2,k 


= 7,11,15... 



5.3 Radial normalization 

Scalar product of two basis vectors for radial oscillators can be defined by a 
generic form of improper volume integral 

= [ + °° fUx)fiAx)e~ x2 dV, (22) 

J 

where dV represents volume of an elementary sphere in N-dimensional space. 
Using general formula for a volume of a sphere with radius x in N dimensions 

where 

r(o) = i 

r(i/2) = v^F 

r(i) = i 

r(jv) = (N — 1)V(N — 1), 

the elementary volume dV, and thus volume integrals fl2~2|) can be easily 
computed. Specificly, the 'volume' of sphere with radius x in dimensions iV = 
1,2,3 are given by 2x, ttx 2 , 4/3ttx 3 , etc. To simplify Haskell computations 
we will ignore constant factors in definition of the N-dimensional volume and 
define dV simply as dV = dx, xdx, x 2 dx, etc., for iV = 1, 2, 3.., respectively. 



5.3.1 One dimensional 

The 'volume' integral in one dimension is twice the value of the one di- 
mensional integral over [0, +oo). This is, generally, not the same as the 
improper integral over (— oo, +oo) - unless the function being integrated is 
even: g(—x) = g(x). 

The recurrence formulae for integrals P n = 2 J +o ° x n e~ x dx, are given 
below. 

Po = 

Pi = 1 

P n - 1 P 

r n — 2 r n-2- 
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As opposed to the case of one dimensional linear normalization the integrals 
P n for odd powers of n do not vanish. 

The case I = admits all solutions k — 0, 1, 2..., but further investigation 
reveals that they do not form the basis since scalar products of one odd and 
one even eigenf unctions is not equal to zero, as in this example, 

(0,0|0,1) = 0.797884. 

However, the set of solutions for I = can be split into two sets: one even 
{k = 0, 2, 4, 6...) and one odd {k = 1, 3, 5...) - each forming the basis 

(0,k\0,k') ^ab- 
solutions for the remaining oscillators (/ > 0) are restricted to odd eigen- 
functions [k = 1,3,5..), because even eigenfunctions are unphysical. As a 
result we can summarize this case as follows: 



Even k basis 



/ = k = 0,2,4... E Qtk = 1,5,9.. . 



Odd k bases 



/ = 





k = 


1,3, 


5... 




= 3,7,11... 


/ = 


1 


k = 


1,3. 


5... 




= 5,9,13... 


/ = 


2 


k = 


1,3. 


5... 


E2,k 


= 7,11,15... 



This table is particularly interesting because it summarizes what is tradi- 
tionally known as three dimensional isotropic harmonic oscillator \TE\. When 
expressed in spherical coordinates and after separation of variables the radial 
part of its Hamiltonian becomes 



2m r dr 2 



1 2 2 l(l + l)tT 

H = —- —r H — mu r + 



2mr s 



which can be further reduced to the form ([!]) after representing the radial 
solutions as 

Ri,k = -<f>i,k( r )- 

r 

Functions 0^ corespond to fi t k£~ r2 ^ 2 investigated in this paper, and the 
three dimensional normalization of solutions R^/. becomes one dimensional 
normalization of (pi k 
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roo roc 

/ Ri, k Ri,k> 4vrr 2 dr = Air / 4>i,k<^i,k' dr. 
Jo Jo 



According to (13) energies are given by Ei^ = 2(l + k) + l, in units of 



Hid 



It is customary to represent this as e n = (n + 3/2)hu, where n = I + k — 1 - 
taking one of the values n = 0, 1, 2, 3... 

With every n, we can associate either one or many pairs of (I, k), which 
generate the states \l, k) coresponding to the same energy e n 



n = =^ |0,1) 

n = l |1,1) 

n = 2 =^ |0,3), |2, 1) 

n = 3 1 1, 3>, |3, 1) 

n = A |0,5), |2,3>, |4, 1) 



The energy levels are degenerate; when n is even then there are | + 1 
states 

|0,n+ 1), |2,n-l), |4,n-3}... 



and when n is odd then there are states 



|l,n), |3,n— 2), |5,n — 4)... 

associated with each energy level. When 21 + 1 possible eigenvalues of L3 are 
also considered for each I then the total degeneracy of energy levels becomes 
|(n+l)(ri+2). The solutions obtained by other means ( see for example [|i~5|1 ) 
are exactly the same as the ones presented in this paper. 



5.3.2 Two dimensional 

According to the condition iV > 21 + 1, we can accept both even and odd 
solutions for oscillator I = 0, while the even solutions for the remaining 
oscillators (I > 0) must be rejected as unphysical. However, it seems that the 
two dimensional normalization does not lead to any well defined pattern of 
orthogonal bases and therefore we reject this ill-specified. But similar 

normalization in the next dimension (N = 3) provides us with several usable 
bases, as shown below. 
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5.3.3 Three dimensional 



In accordance with the general pattern of physical admissibility of even wave- 
functions, N > 21 + 1, the first two oscillators, I = and 1 = 1, admit both 
even and odd solutions k = 0,1,2.... For the remaining oscillators, / > 2, 
only the odd solutions k = 1,3,4... are acceptable. 

For each /, the following pattern of orthogonality relations emerges 

(/, k\l, k') = 5 kk , for \k' - k\ = 0, 4, 8, 12... 

Consequently, the eigenvectors of the first two oscillators can be partitioned 
into four bases - two even and two odd, while each of the remaining oscillators, 
I > 2, has two orthogonal bases. The following four tables summarize the 
three dimensional normalization. 

Even k bases 



I = 

I = 1 


k = 0,4,8... 
k = 0,4,8... 


£ 0)fe = l,9,17... 
£i, fc = 3,ll,19... 




I = 

/ = 1 


k = 2,6, 10... 
k = 2,6,10... 


£ 0l * = 5,13,21... 
£ lifc = 7,15,23... 



The first odd k basis 



I = 





k = 


1,5,9... 




= 3,11,19... 


I = 


1 


k = 


1,5,9... 




= 5,13,21... 


I = 


2 


k = 


1,5,9... 


E2,k 


= 7, 15,23... 



Its degenerated energy levels are given by 

3 

£ n = (n + -)hu, 
where n — I + k — 1 = 0, 1, 2... 

The second odd k basis 



I = 





k = 


3,7,11... 


Eo,k 


= 7,15,23... 


I = 


1 


k = 


3,7,11... 


Ei,k 


= 9,17,25... 


I = 


2 


k = 


3,7,11... 


E 2 ,k 


= 11,19,27... 
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Its degenerated energy levels are given by 

£ n = (n + ^)huj 
where n — I + k — 3 = 0, 1, 2, 3... 

5.3.4 N dimensional 

Specific cases, examined in the last few sections, exhibit certain pattern, 
which can be generalized on N dimensions and directly verified by Haskell 
function bracket. 

In each N dimensional space there exists a treshold integer number I' = 
(N — l)/2, such that for all I < I' all wavefunctions k = 0,1,2... can be 
accepted as physical. Above this threshold the even solutions k = 0, 2, 4 must 
be rejected as unphysical. However, it appears that radial normalization in 
any of the even dimensions N = 2,4,6... does not lead to properly defined 
basis or bases. In contrary, for each of the odd dimensions N = 1,3, 5.. there 
is a clear pattern of several staggered orthogonal bases. Half of them are 
even in k and they are subjected to the limitation of physicality. The other 
half are orthogonal bases with odd vectors k and they are all physical and 
well defined. 

For N = 1 the vectors belonging to the same basis are enumerated by 
Ak = 2. For N = 3 the staggering of bases is given by Ak = 4, for iV = 5 it 
is Ak = 6, etc. 

6 Conclusions 

We have shown that spiked oscillators described by Hamiltonians (|1|) have 
simple and exact eigenfunctions, subject to further normalization restrictions, 
which are specific to a problem dimension. One such restriction is related to 
integrability of the solutions and the other to selection of orthogonal bases. 
The form (|I|) can be interpreted either as a family of one dimensional linear 
spiked oscillators, or families of radial spiked oscillators in odd N dimensional 
spaces, where N = 1,3,5... Even dimensions N = 2,4,6.. must be rejected 
because they do not give rise to orthogonal bases. 

Each Hamiltonian,if/, has an uncountable number of eigenvectors \l,k), 
where k = 0,1,2.. Generally, all solutions indexed by odd k = 1,3,5.. are 
integrable and constitute one basis or several interleaved bases. In contrary, 
only a limited number of oscillators I admit even eigenfunctions k = 0,2,4.. 
for a given odd dimension iV = 1,3,5.. - each forming a basis or a set of 
interleaved bases. 
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Each Hamiltonian (|I]) leads to an equidistant energy spectrum, isomor- 
phic to a spectrum of ordinary quantum oscillator, but subjected to restric- 
tions described above. Our factorization method which is akin to SUSY 
method, and which relies on two kinds of operators: intertwing operators bi\ 
and bi, and ladder operators arf and a h is purely algebraic. But while the 
SUSY "method correlates solutions of two Hamiltonians with isomorphic spec- 
tra, our method extends this approach on infinite set of Hamiltonians. The 
theory is augmented by very simple Haskell program, which directly imple- 
ments these operators, generates the eigenfunctions, tests integrability of the 
solutions and verifies orthogonality conditions for wavefunctions enumerated 
by quantum numbers I and k. 
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